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The  humpback  whale  (. Megaptera  novaeangliae)  is  exceptional  among  the  baleen  whales  in  its 
ability  to  undertake  acrobatic  underwater  maneuvers  to  catch  prey.  In  order  to  execute  these  banking 
and  turning  maneuvers,  humpback  whales  utilize  extremely  mobile  flippers.  The  humpback  whale 
flipper  is  unique  because  of  the  presence  of  large  protuberances  or  tubercles  located  on  the  leading 
edge  which  gives  this  surface  a  scalloped  appearance.  We  show,  through  wind  tunnel  measurements, 
that  the  addition  of  leading-edge  tubercles  to  a  scale  model  of  an  idealized  humpback  whale  flipper 
delays  the  stall  angle  by  approximately  40%,  while  increasing  lift  and  decreasing  drag.  ©  2004 
American  Institute  of  Physics.  [DOI:  10.1063/1.1688341] 


The  turning  performance  of  aquatic  animals  is  con¬ 
strained  by  morphology  with  the  mobility  and  flexibility  of 
the  body  and  the  hydrodynamic  characteristics  and  position 
of  the  control  surfaces  influencing  the  animal’s  turning  per¬ 
formance.  The  baleen  whales  of  the  cetacean  suborder  Mys- 
ticete  are  the  largest  of  all  animals.  Much  of  the  body  is 
inflexible  because  of  their  specialized  feeding  system,  which 
restricts  their  maneuverability.  However,  the  humpback 
whale  ( Megaptera  novaeangliae)  is  unusual  compared  to 
other  mysticetes  in  its  ability  to  undertake  acrobatic  turning 
maneuvers  in  order  to  catch  prey. 1,2  Humpback  whales  uti¬ 
lize  extremely  mobile,  high-aspect-ratio  flippers  for  banking 
and  turning.  Large  rounded  tubercles,  shown  in  the  photo¬ 
graphs  of  Fig.  1,  along  the  leading  edge  of  the  flipper  are 
morphological  structures  that  are  unique  in  this  animal.  The 
position,  size  and  number  of  tubercles  on  the  flipper  suggest 
analogs  with  specialized  leading-edge  control  devices  asso¬ 
ciated  with  improvements  in  hydrodynamic  performance  on 
biological3  and  engineered4  lifting  surfaces.  It  has  been  hy¬ 
pothesized  that  these  leading-edge  tubercles  can  modify  the 
hydrodynamic  characteristics  of  the  flipper  to  increase  its 
effectiveness  in  turning.5'6  In  this  study,  we  report  wind  tun¬ 
nel  measurements  of  the  lift  and  drag,  along  with  other  pa¬ 
rameters,  as  a  function  of  the  angle  of  attack  of  idealized 
scale  models  of  a  humpback  whale  flipper  with  and  without 
tubercles. 

We  constructed  two  scale  models  of  a  humpback  pecto¬ 
ral  flipper;  one  with  and  one  without  leading-edge  tubercles. 
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We  show  both  of  these  models  in  Fig.  2.  The  flipper  plan- 
form  geometry  was  modeled,  roughly,  after  the  left  pectoral 
flipper  of  a  9.02  m  male  humpback  whale.6  Photographs  of 
whale  flippers  reveal  variability  in  the  size  and  shape  of  the 
tubercles  among  individuals.  Therefore,  we  did  not  intend, 
here,  to  model  a  specific  flipper  but  rather  capture  the  gross 
morphology  that  is  common  to  the  species.  We  do  expect  that 
the  details  of  the  geometry  (for  example,  abrupt  changes  in 
surface  curvature)  would  influence  the  performance  of  an 
individual  flipper.  Our  first  (smooth)  model  flipper  had  a 
smooth  leading  edge  with  a  shape  that  closely  matched  the 
average  leading-edge  geometry  of  the  animal  while  our  sec¬ 
ond  (scalloped)  flipper  had  a  sinusoidal  leading-edge  profile 
approximating  the  pattern  found  on  the  specimen  flipper. 
Consistent  with  the  animal,  the  sinusoidal  pattern  had  both 
an  intertubercle  spacing  and  a  tubercle  amplitude  that  de¬ 
creased  with  spanwise  (distal)  location.  The  flipper  cross- 
sectional  profile  is  similar  to  conventional  turbulent-flow  air¬ 
foils  with  a  maximal  thickness  of  approximately  20%. 6 
Therefore,  we  constructed  our  models  based  on  a  symmetri¬ 
cal  NACA  0020  foil  section.  Both  model  flippers,  machined 
on  a  CNC  mill  from  3.81  cm  thick  clear  polycarbonate,  had 
a  maximal  chord  of  16.19  cm  and  a  span  of  56.52  cm.  The 
surface  was  coated  with  self-leveling  epoxy  (West  Systems). 
After  curing,  the  surfaces  were  wet  sanded  with  progres¬ 
sively  finer  emery  paper  down  to  600  grit. 

Our  tests  were  performed  in  the  low  speed  closed-circuit 
wind  tunnel  at  the  United  States  Naval  Academy.  The  tests 
were  conducted  at  incompressible  flow  conditions,  with  the 
aerodynamics  simulating  the  hydrodynamics  of  flow  around 
the  flipper.  As  the  flippers  are  not  oscillated  during  turning 
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FIG.  1.  Left:  photograph  of  a  breaching  humpback  whale  ( Megaptera  no- 
vaeangliae)  showing  the  size  of  the  pectoral  flippers  relative  to  the  animal. 
Right:  detail  view  of  a  pectoral  flipper.  This  flipper  shows  the  characteristic 
tubercles  giving  the  leading  edge  a  scalloped  shape.  The  trailing  edge  is 
relatively  free  from  these  morphological  complexities.  Photographs:  Dennis 
Bowen. 

maneuvers  by  the  whale,  we  were  justified  in  considering 
only  steady  state  flow  conditions  in  testing.  The  maximum 
Mach  number  reached  during  our  tests  was  0.2.  The  vented 
test  section  has  a  cross  section  of  54X38  in.  (137X97  cm) 
and  a  length  of  94  in.  (239  cm).  A  contraction  ratio  of  6.23 


FIG.  2.  Smooth  (left)  and  scalloped  (right)  idealized  humpback  whale  flip¬ 
per  models  used  for  wind  tunnel  testing  of  lift,  drag,  efficiency,  and  stall 
behavior.  While  testing,  only  one  of  the  models  is  mounted  in  the  wind 
tunnel.  These  idealized,  scale  models,  milled  from  clear  polycarbonate 
sheet,  use  a  NACA  0020  sectional  profile  (20%  thick). 


enhances  production  of  a  uniform  flowfield,  as  does  a  qua¬ 
druple  screen  configuration  in  the  upstream  settling  chamber. 
The  wind  tunnel  experiments  were  conducted  at  uncorrected 
Reynolds  numbers  (Re=cf/00r,_1)  in  the  range  5.05 
X  10S-5.20X  105  as  referenced  to  the  freestream  velocity, 
U* ,  mean  chord,  c,  and  the  kinematic  viscosity  of  the  air,  v. 
For  the  whale,  we  would  expect  the  Reynolds  number  to  be 
in  the  neighborhood  of  Re~- 1 06  based  on  a  flipper  mean 
chord  of  50.7  cm6,  a  lunge-feeding  speed  of  2.6  ms-1  and 
seawater  viscosity  of  v=  1.35X  ICE2  cm2  s_1.  Although  the 
Reynolds  number  used  in  our  experiments  is  approximately 
half  of  that  estimated  for  the  animal,  it  is  well  within  the 
operating  range  of  adult  animals  and  perhaps  close  to  the 
maximum  for  young  animals.  Nevertheless,  we  performed 
tests  on  the  effect  of  the  Reynolds  number  on  the  measured 
quantities  and  found  that  for  Re>4X105,  the  lift  coefficient 
was  relatively  insensitive  to  Re  at  moderate  incidence  angles. 

Each  flipper  model  was  mounted  vertically  in  the  wind 
tunnel  on  a  rotating  yaw  table  which  changed  the  local  pitch 
angle.  The  models  were  rotated  through  an  angle  of  attack, 
a,  range  of  —2°  to  +20°  (A a  =  0.5°)  using  a  pitch-pause 
technique.  This  included  some  hysteresis  tests  that  involved 
first  increasing  the  angle  of  attack  from  —2°  to  +20°,  and 
then  decreasing  the  angle  of  attack  from  +20°  to  —2°. 

Aerodynamic  forces  and  moments  were  transmitted  to  a 
six-component  Schenck  external  compact  platform  balance 
located  beneath  the  test  section  yaw  table.  The  standard  de¬ 
viation  in  force  measurement  over  a  representative  15- 
minute  test  at  a  constant  flow  condition  has  been  demon¬ 
strated  to  be  ±0.03  lb  (0.13  N)  in  drag,  ±0.24  lb  (1.08  N)  in 
lift,  and  ±  0.08  ft-lb  (0. 1 1  Nm)  in  pitching  moment.  The  co¬ 
efficients  for  these  forces  have  standard  deviations  of 
±0.0003,  ±0.003,  and  ±0.001,  respectively,  for  the  present 
wind  tunnel  models.  Additionally,  the  test  section  was  instru¬ 
mented  with  pressure  transducers  and  a  thermocouple  to  ob¬ 
tain  real-time  measurement  of  freestream  velocity  and  air 
density. 

The  millivolt-level  balance  output  signals  were  moni¬ 
tored  and  acquired  using  a  6y-digit  Agilent  34970 A  data 
acquisition/switch  unit  with  a  34901A  multiplexer  module. 
Serial  connections  to  a  Jetter  Mikro  PLC  were  accessed  for 
monitoring  the  model  orientation  encoders.  A  conventional 
HP-IB  (IEEE  488.2)  communications  link  interfaced  the  data 
acquisition  unit  with  the  host  laptop  computer.  Scan  rates  as 
high  as  250  channels  per  second  were  executed  with  this 
system. 

During  data  acquisition,  nine  channels  were  scanned  (six 
balance  outputs  including  three  force  measurements  and 
three  moment  measurements,  the  balance  excitation  voltage, 
the  test  section  dynamic  pressure,  and  the  freestream  air  tem¬ 
perature).  A  single  scan  was  taken  at  6y-digit  precision  over 
an  integration  aperture  of  20  power  line  cycles,  as  prescribed 
for  high  common  mode  rejection  on  low-level  signals.  The 
effective  scan  rate  was  0.6  channels/s.  Approximately  5  s 
were  allocated  for  tunnel  transients  to  diminish  prior  to  trig¬ 
gering  the  data  acquisition  system.  Each  data  acquisition 
scan  process  required  about  6  s  to  complete.  Further  details 
of  the  measurement  system  are  available  elsewhere.7 
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FIG.  3.  Lift  and  drag  data  for  the  smooth  whale  flipper  model,  (a)  Wind 
tunnel  measurement  of  the  lift  coefficient,  CL ,  shown  as  a  function  of  the 
angle  of  attack,  a.  (b)  Drag  coefficient,  C D ,  shown  as  a  function  of  the 
angle  of  attack,  (c)  Aerodynamic  efficiency,  LID ,  which  quantifies  the  drag 
cost  of  producing  lift,  shown  as  a  function  of  the  angle  of  attack,  (d)  Plan- 
form  profile  of  the  smooth  flipper  model  (solid  line)  superimposed  on  the 
actual  flipper  profile  (dotted  line).  In  this  orientation,  flow  would  progress 
from  top  to  bottom,  (a)  and  (b)  show  that  stall  occurs  near  a***  ±12°. 


The  primary  nondimensional  parameters  describing  the 
aerodynamic  forces  placed  on  the  models  by  the  flow  are  the 
drag  coefficient,  CD=2D(pUl0A)~l,  and  the  lift  coefficient, 
CL=2L{pU^A)  .  While  we  also  measured  the  side  force 
and  the  three  moment  coefficients,  we  do  not  report  these 
here.  In  the  expressions  for  the  drag  and  lift  coefficients,  D 
and  L  are,  respectively,  the  drag  force  and  the  lift  force,  and 
p  is  the  air  density.  The  values  of  the  projected  planform 
area.  A,  and  the  mean  chord,  c,  for  both  the  smooth  and 
scalloped  models  are  737.7  cm2  and  13.04  cm. 

Reynolds  number  sensitivity  was  established  using  the 
smooth  flipper  model.  With  the  Reynolds  number  varied 
from  135  000  to  550  000,  the  lift  and  drag  coefficients  were 
measured  at  angles  of  attack  of  0°  and  5°.  Although  both  the 
lift  and  drag  coefficients  show  a  notable  dependence  at  the 
lower  Reynolds  numbers,  the  curves  flatten  out  to  indicate 
very  little  or  no  Reynolds  number  dependence  at  the  test 
conditions  of  the  present  study.  The  lift  coefficient  is  the 
most  insensitive,  with  essentially  no  dependence  for  Rey¬ 
nolds  numbers  above  200000.  The  drag  coefficient  ap¬ 
proaches  a  constant  value  for  our  test  angles  of  attack  above 
Reynolds  numbers  of  approximately  400000,  although  a 
truly  constant  trend  is  not  established.  We  expected  that,  be¬ 
cause  of  the  narrow  tip,  the  stall  behavior  on  the  outboard 
portion  of  the  flipper  would  be  strongly  dependent  on  the 
Reynolds  number,  especially  at  the  lower  range.  Neverthe¬ 
less,  the  CD  values  decrease  at  such  a  low  rate  at  the  upper 
range  of  the  Reynolds  number  sweep,  particularly  for  a 
=  5°,  that  the  Reynolds  sensitivity  can  be  assumed  to  be 
negligibly  small. 

The  lift  curve  for  the  smooth  model.  Fig.  3(a),  shows  a 
linear  relation  between  the  lift  coefficient,  CL ,  and  the  angle 
of  attack,  a,  for  angles  below  12°.  In  this  region,  the  slope  of 
the  lift  curve  has  the  value  a=dCL/da  =  Q.0169  deg-1. 
This  measured  slope  is  within  5%  of  the  results  calculated 


from  thin  airfoil  theory8  (based  on  an  effective  aspect  ratio  of 
8.7).  The  constant-slope  segment  of  the  lift  curve  indicates 
that  the  flow  is  attached  to  the  surface  of  the  model  flipper 
and  stall  is  absent  during  this  flow  state. 

The  onset  of  stall  occurs  suddenly  but  the  loss  in  lift  is 
not  as  considerable  as  we  would  expect  with  complete 
leading-edge  flow  separation.  For  the  smooth  flipper  model, 
“stall  occurs  at  12°,  yet  the  drop  in  CL  is  only  13%.  This 
small  loss  in  lift  might  be  due  to  flow  separation  near  the 
distal  region  before  flow  separation  near  the  proximal  region 
resulting  in  only  a  partial  loss  of  lift.  While  we  did  not  per¬ 
form  flow  visualization  to  validate  this  hypothesis,  the  high 
taper  of  the  outboard  portion  of  the  flipper  suggests  that  low 
Reynolds  number  effects  are  present  here  that  might  not  be 
present  at  inboard  stations.  These  low  Reynolds  number  ef¬ 
fects  could  result  in  partial  stall  at  the  distal  end  of  the  flip¬ 
per.  The  subsequent  decrease  in  CL  from  0.78  to  0.38  reveals 
a  progression  of  stall  as  we  increase  a  from  12.1°  to  18.5°. 

The  drag  coefficient,  CD  [Fig.  3(b)],  shows  that  the  onset 
of  stall  occurs  with  an  increase  in  drag  from  CD  =  0.0485  to 
CD  =  0.0821.  The  69%  increase  occurs  over  a  change  in  a  of 
approximately  0.5°.  Then  CD  increases  more  gradually  from 
0.0821  to  0.2226.  This  2.7-fold  increase  in  the  drag  indicates 
a  progression  to  total  flow  separation  over  the  entire  flipper 
at  greater  a. 

The  LID  ratio,  or  aerodynamic  efficiency,  shows  a  broad 
range  of  operation  at  a  high  LID  as  shown  in  Fig.  3(c).  This 
curve  characterizes  the  drag  cost  of  producing  lift.  The  peak 
occurs  at  a  =  7.5°  with  a  value  of  (L/Z))max= 22.5.  This  is 
rather  efficient  for  a  section  having  a  20%  thickness.9  More 
significantly,  the  smooth  flipper  maintains  LID  values  over 
20  through  a  broad  range  of  incidence  angles:  a  =  5-10°. 

We  now  shift  our  focus  to  the  aerodynamic  performance 
measures  of  the  scalloped  flipper  model  that  show  that  the 
presence  of  the  tubercles  delays  stall  to  higher  a,  increases 
maximum  lift,  and  decreases  drag.  For  example,  below  an 
incidence  of  approximately  8.5°,  the  lift  curve  is  largely  un¬ 
changed  from  the  smooth  flipper  model  results  [Fig.  4(a)]. 
However,  the  lift  curve  takes  a  distinctly  shallower  slope  of 
a  =  0.494  deg-1  from  a  =  8.5°  - 14.5°,  yet  the  trend  remains 
largely  linear.  Although  the  CL  for  the  scalloped  model  is 
lower  than  the  CL  for  the  smooth  flipper  for  a  limited  range 
(9.3^a,«12),  it  is  equal  to  or  greater  than  the  results  from 
the  smooth  model  at  other  angles.  We  also  observe  that  «stan 
increases  by  40%  to  16.3°.  Of  significant  interest  and  utility 
is  the  fact  that  the  value  of  C,  increases  6%  over  the 

max 

smooth  model,  from  0.88  to  0.93.  Thus,  the  lift  of  the  scal¬ 
loped  flipper  model  shows  an  expanded  operating  envelope: 
higher  lift  at  higher  incidence  angles. 

The  drag  curve  for  the  scalloped  flipper,  shown  in  Fig. 
4(b),  shows  pronounced  increases  with  the  onset  of  each 
aerodynamic  state  (for  example,  attached  flow,  partial  stall, 
and  deep  stall).  At  low  angles,  the  drag  is  no  greater  than  the 
smooth  results.  The  first  break  in  the  lift  curve  is  accompa¬ 
nied  by  a  drag  increase  of  A  C D  =  0.0 1 07  or  34%,  making  the 
drag  higher  than  the  smooth  results  for  the  limited  range  of 
10.3<  a<  1 1.8.  Beyond  this  range,  however,  the  scalloped 
flipper  model  produces  consistently  lower  drag  than  the 
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FIG.  4.  Lift  and  drag  data  for  the  scalloped  whale  flipper  model.  The  solid 
lines  in  (a),  (b),  and  (c)  show,  for  comparison,  the  average  of  the  data  from 
the  wind  tunnel  experiments  utilizing  the  smooth  flipper  model  [Figs.  3(a)- 
3(c)],  (a)  Wind  tunnel  measurement  of  the  lift  coefficient,  CL ,  shown  as  a 
function  of  the  angle  of  attack,  a.  Note  that  the  scalloped  geometry  main¬ 
tains  lift  to  higher  angles  of  attack  than  the  smooth  model,  (b)  Drag  coeffi¬ 
cient,  CD  ,  shown  as  a  function  of  the  angle  of  attack,  (c)  Aerodynamic 
efficiency,  L/D.  (d)  Profile  of  the  idealized  scalloped  model  (solid  line) 
superimposed  on  the  actual  whale  flipper  profile  (dotted  line).  These  mea¬ 
surements  show  that  the  scalloped  model  delays  complete  stall  to  asta]1 
±  17.5°,  a  40%  improvement  over  the  stall  angle  of  the  smooth  model. 


smooth  model,  by  as  much  as  32%.  Therefore,  the  expanded 
operating  envelope  is  accompanied  by  a  significant  reduction 
in  drag,  meaning  the  animal  can  expend  less  energy  in  gen¬ 
erating  greater  maneuvering  forces  with  the  flipper. 

As  a  result  of  the  lift  and  drag  characteristics,  aerody¬ 
namic  efficiency  is  decremented  at  some  incidence  angles, 
but  incremented  at  most  others.  Figure  4(c)  shows  LID  for 
the  scalloped  geometry.  The  peak  aerodynamic  efficiency  oc¬ 
curs  around  the  same  point  as  the  smooth  case  {a  =  1.5°)  but 
has  a  higher  value  of  23.2.  However,  the  range  of  highest 
values  is  reduced  somewhat:  the  scalloped  model  has  a  lower 
LID  for  a  =  5°-8.5°.  The  LID  ratio  is  greater  at  all  other 
incidence  angles.  This  means  that  the  scalloped  flipper  per¬ 
forms  better,  proportionately,  at  most  points  in  the  opera¬ 
tional  envelope,  particularly  at  higher  angles. 

Our  findings,  specifically,  are  that  the  scalloped  leading 
edge  of  the  flipper  serves  to  delay  stall  by  providing  higher 
lift  at  higher  incidence  angles,  and  to  ameliorate  the  post¬ 
stall  characteristics  by  maintaining  higher  lift  with  lower 
drag.  The  ability  to  augment  performance  at  high  angles  of 
attack  would  be  advantageous  for  the  humpback  whale  while 
feeding  and  maneuvering  where  an  expanded  operating  en¬ 
velope  is  needed.  For  example,  the  imbalance  in  lift  between 
the  flippers  (high  positive  lift  on  one,  high  negative  lift  on 
the  other)  creates  an  enhanced  ability  to  roll,  and  to  do  so  at 


higher  rates.10  Maximum-lift  orientation  caused  by  placing 
both  flippers  in  the  same  orientation  creates  a  large  pitching 
moment  about  the  animal’s  center  of  mass,  enabling  the  ani¬ 
mal  to  dive  or  climb  quickly.  Additionally,  a  greater  imbal¬ 
ance  in  drag  gives  the  whale  added  yaw  control  that  allows 
the  whale  to  turn  swiftly.  Enhanced  maneuvering  capabilities 
permit  tighter  turns,  which  give  these  whales  an  advantage  in 
their  use  of  “bubble  netting”  to  catch  prey.1  2,6  We  should 
also  note  that  the  animal  might  use  “controlled”  stall  to  en¬ 
hance  maneuvering. 

From  the  hydrodynamic  performance  capabilities  of  the 
flippers  with  the  addition  of  the  leading-edge  tubercles,  we 
can  draw  analogies  to  vortex  generators  that  energize  the 
flow  over  an  aircraft  wing.  These  vortex  generators,  usually 
small  tabs  placed  near  the  leading  one-third  of  the  chord  of 
an  aircraft  wing,  cause  a  greater  momentum  exchange  within 
the  boundary  layer  than  would  occur  without  the  vortex  gen¬ 
erators.  This  additional  momentum  transport  helps  to  keep 
the  flow  attached  to  the  lifting  surface  despite  an  adverse 
pressure  gradient.  The  dynamics  of  the  vortices  serve  to 
maintain  or  increase  lift  by  preventing  stall.4  The  lift-force 
and  drag-force  modification  created  by  a  scalloped  leading 
edge  produce  results  similar  to  these  produced  by  vortex 
generators.  Therefore,  the  enhanced  lift  and  drag  effects  af¬ 
forded  by  a  scalloped  leading  edge  offer  a  possible  alterna¬ 
tive  to  tab-like  vortex  generators  for  specialized,  high- 
performance  lifting  surfaces  for  low  Reynolds  number 
autonomous  aircraft  and  for  underwater  vehicle  operation. 
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